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We present a novel procedure for calculating non-equilibrium two-point Green's functions in the 
O(N) cj) 4 theory at large N. The non-equilibrium density matrix p is constructed via the Jaynes- 
Gibbs principle of maximal entropy and it is directly implemented into the Dyson- Schwinger equa- 
tions (DSE) through initial value conditions. In the large N limit we perform an explicit evaluation 
of two-point Green's functions for two illustrative choices of p. 
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I. INTRODUCTION 

In the present work we show a particular approach to 
the non-equilibrium QFT dynamics based in the Jaynes- 
Gibbs principle of maximal entropy Wa . In contrast to 
other methods in use ua the Jaynes-Gibbs method con- 
structs a density matrix p directly from the observed 
macroscopic quantities (e.g. pressure, density of en- 
ergy, magnetisation, particle current, local momentum, 
local angular velocity, ionisation rate (if plasma is in 
question), etc.). The p is then implemented through 
the generalised Kubo-Martin-Schwinger (KMS) condi- 
tions into the dynamical equations for Green's functions. 
To keep complexity minimal we illustrate our method on 
a paradigmatic physical system described by A0 4 theory 
with O(N) internal symmetry. The plan of this paper is 
as follows. 

In Sec. II we review the Jastpes-Gibbs principle. The 
role of the Shannon entropy TO is emphasised. Sec. Ill is 
devoted to the construction of the DSE for QFT systems 
away from equilibrium. We use the canonical formalism 
which appears to be very natural in this context. In order 
to reflect the density matrix in the dynamical equations 
we show how to formulate the relevant initial-time condi- 
tions. The DSE for two-point Green functions are worked 
out in Sec. IV. With this mathematical setting, we take 
in Sec.V the large N limit. In the latter case the DSE 
for two-point Green functions are decoupled (virtually 
by chance). We explicitly solve these for two illustrative 
choices of the initial-time conditions. 



II. INITIAL CONDITIONS, JAYNES-GIBBS 
PRINCIPLE OF MAXIMAL ENTROPY 

In this section we would like to briefly review the 
Jaynes-Gibbs principle xiLmaximal entropy (also maxi- 
mum calibre principle) tn3ll2l . The objective of the prin- 
ciple is to construct the 'most probable' density ma- 
trix which fulfils the constraints imposed by experimen- 
tal/theoretical data. 

The standard rules of statistical physics allows us to 
define the expectation value via the density matrix p as 



(•••) = Tr(p • • •), with the trace running over a complete 
set of physical states describing the ensemble in question 
at some initial time ti. _ _ 

The usual approaches Wu trying to determine p 
start with the Liouville equation and hence with the 
Schrodinger picture. Rather than following this path, 
we shall use the Heisenberg picture instead. This will 
prove useful in Sec. III. 

In order to determine p explicitly we shall KWpt to the 
Jaynes-Gibbs principle of maximal entropy HErEiJ . The 
basic idea is 'borrowed' from the information theory. Let 
us assume that we have criterion of how to characterise 
the informative content of p. The most "probable" p is 
then selected out of those p which are consistent with 
'whatever' we know about the system and which have 
the last informative content (Laplaces's principle of in- 
sufficient reasons). 

It remains to characterise the information content 
(measure) I[p] of p. This was done by C.Shannon tf with 
the result: I[p] = Tr(p hip). 

The density matrix^] is then chosen to minimise I[p\. 
Note that no assumption about the nature of p was made; 
namely there was no assumption whether p describes 
equilibrium or non-equilibrium situation. To_jpiit more 
flesh on the bones, let us rephrase the former tin. What 
we actually need to do is to maximise Sg subject to the 
constraints imposed by our knowledge of expectation val- 
ues of certain operators Pi [<$, d(p] , . . . , P n [(f>, dej)] . In con- 
trast to equilibrium, all Pk [. . .] 's need not to be the con- 
stants of the motion. So namely if one knows that 



(P k [(p,d4>]) = f k (x 1 ,x 2 , . • .), 
the entropy maximalisation leads to 

P= Z[Ai,...,A„] 



(1) 



* It can be shown H that —I[p] (also called the informa- 
tive entropy) equals (in base 2 of logarithm) to the expected 
number of binary (yes/no) questions whose answer takes us 
from our current state of knowledge to the one certainty. 



1 



with the 'partition function' 

Z[X U . . . , A„] = Tr (, E HI *-> *<n,...W*M) 

It is possible to show that the stationary solution of So is 
unique and maximal H. In the previous cases the time 
integration is not either present at all (fk is specified 
only in the initial time ti), or is taken over the gathering 
interval (— r, ti). 

The Lagrange multipliers A& might be eliminated if one 
solves n simultaneous equations 



SlnZ 



SX k (xi, ...)' 

The explicit solution of ^ may be formally written as 

S Sc[fl, ■ ■ ■ , fn]\ max 



(2) 



AfcOi, . . .) 



Sfk(xi,...) 



III. OFF-EQUILIBRIUM DYNAMICAL 
EQUATIONS 

In this section we derive the off-equilibrium Dyson- 
Schwinger equations using the canonical formalism. We 
believe that this is a new and far more natural formu- 
lation for systems away from equilibrium. The more i 
tuitive path-integral formulation of Calzetta and Hu t 
is not particularly suitable-in this case, because the con- 
nection between kernels 13 and initial-time constraints 
turns out to be rather non-trivial O. 

Let us deal first with a single field <f>. We start with 
the action S where <f> is linearly coupled to an external 
source J(x). For the fields in the Heisenberg picture, the 
operator equation of motion can be written as 



SS 
5ct>{x) 1 



= c/> J } + J(x) = 0, 



(4) 



where the index J emphasises that cj> is implicitly J- 
dependent. It will prove useful in the following to re- 
express <f> in such a way that the J-dependence will 
become explicit. The latter can be done via an unitary 
transformation connecting <fi J (governed by H — Jcf>) with 
4> (governed by H). If J(x) is switched on at time t = ti 
we have 



/, 



5 SI 



J 



exp(z I d y J{y)4>{y)) 
'c 



0, (5) 



with T* being the T*-ordering. The close-time path C is 
the standard Keldysh-Schwinger path. Associating with 
the upper branch of C index '+' and with the lower one 
the index ' — ' one may introducing the metric ((13)0/3 ((73 
is the usual Pauli matrix and a,j3= {+; — }) and write 
- J_0_ = J a {a 3 ) a P <f,p = J a (a 3 ) a (3 <fP. For 



the raised and lowered indices we simply read: + = + 
and <p- = —(f>~ (similarly for J a ). Taking Tr(p. . .) with 
P = p[<P,9(f), . . .], we get 
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SS 



Z[J] 5<j>(x) 



SJ a (x) 



Z[J] = -J a (x), (6) 



with Z\J\ being the generating functional for Green's 
functions. Because of the T*-ordering the time deriva- 
tives could be pulled out of [...]. Eq.(^J) may equivalcntly 
be written as 



J a (x) 



SS 

8<j){x) 



SJ a (x) 



(7) 



The 1 indicates the unit. Analogously as for equilibrium 
systems, we have defined the mean field <f> c a (x) as 



(3) cf) ca {x) = {(j) a [x))j = (o- 3 ) Q| a 



5W[J] 
5Jp{x) 



Z[J] = exp(iW[J]). 



Summation over contracted indices is understood. The 
effective action T[cj) c ] is connected with W[J] via the Leg- 
endre transform: T[cf) c ] = W[J] — J c d 4 y J(y)(j) c (y). 

With this mathematicaLseiting we obtain the usual 
equilibrium-like identities oH 



sr\ 



6<f> ca (x) 



= -J a (x), 



d A y G a p{ Xl y) (a 3 Y 3S T^(y,z) = (a 3 ) aj S 4 (x - z), 



with 



-G a p{x,y) = i{T c {4> a {,x)(f) fj {y)}) - i{cf> a (x))((f) (y)) , 



S 2 T 



r$(z,y)- 



For the physical process (i.e. J± — 0) we have from (0) 

<sr[0 c ] ss 



t>?(x) 5<f>(x) 



{4>ca(x) 



+ i / d*yG afi (x,y) (a 3 )^ 



1=0. (8) 



It is worthy of noticing that the LHS of (||) offers a direct 
prescription for a calculation of ST[(j) c ]/S(p"(x). 

So far we have not taken into account the constraints. 
This can be done quite simply. One just sets in p to 
be the solution of Eq. (||) . Using the identity 

00 

e A Be~ A = y2-C n , C = B,C n = [A,C n ^}, 
£ — ' nl 

n=0 

we get then the generalised KMS condition 
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oo 1 

G+-{x, y) = G.+foy) -rTty, 4>{x)C n {y)). 

n=l 

Here A = hip, B = <p(y) and yo =ti- Similarly we could 
derive the generalised KMS conditions for the higher 
point Green's functions 

Eq.(^|) and its successive J variation provide us with 
the coupled integro-differential equations. As a result, 
we get an infinite hierarchy of coupled equations which, if 
furnished with the corresponding KMS conditions, con- 
stitute, in principle, a complete description of the be- 
haviour of a non-equilibrium system. 



IV. THE 0{N) <j> 4 THEORY 

Let us illustrate the aforementioned formalism on the 
O(N) <j) A theory. 

The O(N) 4> 4 theory is described by the bare Lagrange 
function 



i N \ ( N V 

^ = ^E(w Q ) 2 --o(^) 2 )-^ E(^ Q ) 2 

a=l \a=l ) 

Using the explicit form (^|), Eq. (||) reads 
ST 



(9) 



b aa (x) 



'2N 



N 



N 



c («) (€ W) 2 + i<& & E G - x ) 



.6=1 



6=1 



N 



+i 



6=1 



+ f dS d±w d*z J2 GL 3 l: bb (y, w,z)\ = 0. 
J b=l J 

A successive variation with respect to J{y) generates the 
DSE for the two-point Green's functions. 

The dynamical equations can be significantly sim- 
plified provided that both the density matrix and the 
Hamiltonian are invariant against rotation in the N- 
dimensional vector space of fields. This fact leads 
straightforwardly to the following Ward's identities 111 



SW[J] 



8JZ(z) 



(z) = 0, Vflt, 



5 2 W[J] 



8J aa (x)SJP b (y) 

^(3) ab 



J=0 



.7=0 



Gf ff (x,y) = 8 ab G a p(x,y), 



r a/3 7 (yij2/2,2/3) = 0, Va,6,c, 



al^s^iy^y^ ys, yt) <x 8 ab 8 cd + 8 ac s bd + s ad 8 bc . 



With these results we may write the evolution equation 
for the two-point Green's function as follows 



2 , A . N + 2 



D + m a + y i N G aa (x,x) ) G a p(x,y) 



Xo N + 2 



2 N 



( G< ^l a0 (x,x,x,yfj = -8(y-x)(a 3 ) a0 , (10) 



In the following we shall confine ourselves only to such 
situations where both p and H are O(N) invariant. 



V. G aB {X,Y): O(N) <p 4 THEORY IN THE LARGE N 
LIMIT 



One can show BOO by a detailed study of the large-N 
approximation (virtually using only the Ward's identities 
and properties of T an W) that the vertex functions I^ 2 ™) 
must be of order TV 1- ™. The latter suggests that in the 



dynamical equation (10) we can neglect T^ 4 ' terms. 

Let us mention one more point. If we perform the 
expectation value of the Lagrange function (^|) we find 
that this does not depend on G^ in the TV — ► oo limit, 
indeed 



Tr(p£(x)) = -iN {d ax d°G{x,y)\ x=y -m 2 G(x,x)} 



i\ (N+l) 



{G (4) (x,x,x,x) + {G{x,x)) 2 }. 



Here we have used the fact that Galaa — 3G^ bb = 3G^ . 
So we may equally start with the Lagrange function 

JV 

2 

Ao 



1 N 

^-^((afWf-m^fW) 2 ) 



a=l 



■Y,(r(x)) 2 G(x,x). 



(11) 



a=l 



The former fulfils the identity (£) = (£}|jv^oo. It is also 
simple to see that the DSE derived from £ reads 



D x + m, 



o 



G aa {x,x) ) G a p(y,x) 



-8(x - y)(<T 3 ) a p- 



(12) 



This is precisely the same which one would obtain if the 
large N limit would be performed in the original DSE 

Let us now show how to compute G a p(x,y) for some 
familiar choices of the initial-time constraints. 

(i) Equilibrium 

In this case the constraints are usually chosen to be 
the integrals of motion. The only available integral of the 
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motion is the full Hamiltonian H, and the corresponding 
constraint reads oEf 

{P k [<t>,d<t>))\ ti = (H) = f dT'Cv(T') = F(T), (13) 
Jo 

where Cy is the heat capacity at constant volume V. The 
density matrix p maximising Sg is the density matrix of 
the canonical ensemble: p — ^^zJW~' ^ ne Lagrange 
multiplier (3 is determined from Eq.(p|): 

' dF{T) \dT) v \ dT ) v T 1 ' 

In this case the KMS condition is the well known relation 

G+_(x,t 4 ;y,0) = G_+(x, i 4 - z/3; y, 0). (15) 

The DSE are those in (|l|) with G(x, y) = G(x - y). The 
solutions are the equilibrium propagators in the Keldysh- 
Schwinger formalism, i.e. 

iG ±± {k) = ±l 2 - . + 2 7 r/(|A ;o |) < 5(/c 2 - m 2 r ) 
Ar — mf. ± is 

iG ±T (p) = 2n {6(Th) + f(\k \)} S(k 2 - m 2 r ) , (16) 

where fix) = (exp(f3x) — 1) _1 is the Bose distribution, 
and m 2 = m 2 , + ^-G(0) is the renormalised mass. 

(ii) Non- equilibrium: translationally invariant G a p 

The DSE in the former example were immensely sim- 
plified due to the translational invariance of Green's func- 
tions. If we retain the translational invariance this sim- 
plicity will be preserved also to non-equilibrium. As 
an example, let us choose the following initial-time con- 
straint: 

(P[<f>, dcj>])\ ti = (W(k)> = (W(k)) = , 9 (k), (17) 

where TL means the effective Hamiltonian density derived 
from C The density matrix reads 

_ cxp(-/ d 3 k(3(k)H(]<i)) _ exp(- J d 3 k(3(k)ala k ) 

with /3(k) = (3(k)2^/k 2 + m 2 r and 0Qt) = ^gy. The 
former indicates that the different modes have different 
'temperatures'. The DSE agree with those in the equi- 
librium case. The generalised KMS conditions are 

G+_(fc) = e-^ k °/ 2 G- + (k) (18) 

The solution of ( |l2| ) furnished with (|l8|) formally coin- 
cides with the solution (|l6|). The only proviso is that 

f(\ko\) el 3(k)/2l • 



VI. CONCLUSIONS AND OUTLOOK 

We considered the Jaynes-Gibbs principle of maxi- 
mal entropy. This allowed us to construct the non- 
equilibrium DSE. For the 0(N) 4> 4 theory in the N — » oo 
limit we have explicitly evaluated propagators for two 
choices of translationally invariant density matrices. A 
notable advantage of this approach is that one can get 
the DSE without gome, through the Cornwall- Jackiw- 
Tomboulis formalism 1~Ll which would be formidably dif- 
ficult particularly for more than one constraint. 

The O(N) 4> A theory in the large N limit is a nice toy 
model allowing in many cases to approach the dynamical 
equations analytically. For suitable choices of the trans- 
lationally non-invariant initial-time constraints the sim- 
plicity of the equations is such that one may solve thcsa 
exactly. A more detailed report will appear elsewhere . 
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